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Abstract 

New phenomenon of temporal oscillations of nonlinear Faraday rotation in a 
driven four-level system is predicted. We show that in this system with one 
upper level, under the conditions of electromagnetically induced transparency 
created by a strong coupling field, the polarization rotation of weak probe 
light exhibits slowly damped oscillations with a frequency proportional to the 
strength of an applied magnètic field. This opens up an alternate way to 
sensitive magnetometric measurements. Applications in low-light nonlinear 
òptics such as photon entanglement are feasible. 
PACS numbers: 42.50.Gy, 07.55.Ge 



In the last years there has been growing interest in the study of nonlinear Faraday ro- 
tation (NFR) m in coherent atòmic media. Both the regime of weak fields, low atòmic 
coherence and low optical density |||| and the regime of strong field, large long-lived coher- 
ence and optically dense medium |4j] have been extensively investigated. If in the first case 
the ultranarrow magneto-optic resonances were observed with relaxation rate of 27rxlHz 
||, in high light-power regime the large atòmic coherence giving rise to electromagnetically 
induced transparency (EIT) || allows to enhance essentially the nonlinear Faraday signal in 
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optically thick medium. These experiments J||§, as well as the theoretical estimations J7j, 
demonstrate the significant potential of both techniques to prové substantially the sensitivity 
of NFR-based magnetometers to small magnètic fields. 

An important feature of coherent atòmic medium is its extremely high dispersion 
near sharp atòmic resonance lines such that a large NFR occurs even at very weak magnètic 
fields. The simplest system displaying this feature is a three-level A-atom with ground- 
state Zeeman sublevels interacting with the right (RCP) and left (LCP) circularly polarized 
components of optical light in the presence of a static magnètic field. In usual terms of 
induced-coherence treatment, the RCP and LCP components of the light can be repre- 
sented as playing the role of coupling and probe fields with respect to each other. Having 
equal intensities they excite a maximal Zeeman coherence in high-power regime. On the 
other hand, the steep dispersion at probe frequency in more general case of strong coupling 
and weak probe fields has been used earlier jT(J to detect small magnètic fields by 



measuring the probe phase shift via optical interferometry. These investigat ions were aimed 
at improvement of sensitivity to magnètic fields compared to conventional optical pumping 



magnetometers [I 



We study here, for the first time, the magneto-optical rotation in a strongly driven 
coherent medium. The model we are considering is a four-level system (Fig.l) where, as 
before, the RCP and LCP components of weak probe light E p of frequency u p , tuned to the 
|1>— » |3>and |2>— * |3>transitions respectively, create a Zeeman coherence, while a strong 
coupling field E c with frequency uo c drives the atom on the auxiliary transition |0>^ |3>. 
It is obvious, that in Faraday geometry the Zeeman sublevels |l,2>and the state |0>belong 
to different hyperfine components of the ground state. Moreover, the coupling field E c 
should be circularly polarized. Our motivation for the present work is related to the novel 
features of Faraday rotation arising in this system in the EIT-regime. We show here that the 
Faraday rotation angle oscillates in time with a slowly damped amplitude and a frequency 
proportional to strength of magnètic field. For weak probe light and small magnètic fields 
the damping rate of the oscillations is determined by the lifetime of Zeeman coherence, 



which can be made very long compared to the decay time of excited state . This opens up 
a promising possibility to detect the small magnètic fields via the measurements of Faraday 
signal oscillations with a sensitivity not less than already reached one. 

The underlying mechanism of the long-lived oscillations of Faraday signal is an 
existence of stable coherent superposition of two ground-state sublevels |l>and |2>with 
different energies = =F<5,where 8=giiB^>/h is the Zeeman level shift induced by a magnètic 
field B (g is a Lande factor). In the absence of magnètic field the strong coupling field drives 
the àtoms into population-trapped (or noncoupled) states — çf-\0))~ i— 1,2 (fi c and 
Q p are the Rabi freqüències of the coupling and probe fields: Q c » Q p ), which are mixed 
with each other due to two-photon interaction with the circular components of the probe 
light. This occurs at the time scale of lifetime of upper state |3>. For all time thereafter 
the atom is in a steady-state and does not interact with the fields. When a weak magnètic 
field is applied, the system evolves in the same scenario with except that now the states 
|l>and |2>acquire different phases due to their different energies. This means that in 
the EIT-regime the coherent superposition of the states |1) exp(— iSit) and |2) exp(— i^t) 
is eventually created with the constant probabilities except for a very slow decay due to 
Zeeman coherence damping. These superposition states exhibit the quantum beats in time 
evolution of the system, which are revealed in NFR in the form of slowly damped temporal 
oscillations. Thus, the key point is that to a good approximation the time behavior of 
Zeeman coherence is the oscillatory free evolution, but at the same time the susceptibility 
of probe transition, although small, is different from zero that makes it possible to detect 
these oscillations in Faraday signal. 

We now turn to the analytical calculations of this effect in the scheme depicted in 
Fig.l. The total interaction-picture Hamiltonian corresponding to this system has the form 

H = ftAi|l)(l| + frA 2 |2>(2| + ftA c |3)(3| - fi(íï c |3)(0| + O l· |3> <1| + 2 |3> <2| + H..c.) (1) 

where the Rabi freqüències of coupling field and polarization components of probe light are 
defined as Q c = nzoE c /h and Qi,2 = Vl, p exp{i^i^), Çl p = fiE p /h, ^ is the dipole moment of 
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the i-k transition, fj, = /i 3 i = /i 32 , the amplitudes E CjP are considered real and $i )2 are the 
phases of probe polarization components. The detunings 

A c = lo 30 -u c - S , Ai = lü 3í - Up - Si, i = 1, 2, (2) 

include the frequency shifts of atòmic transitions S lj2 and 5 = goUsBrno/h caused by a 
magnètic field, where go and m are the giromagnetic factor and magnètic quantum number 
(momentum projection) of the level |0). The thermal motion of the àtoms with velocity v 
in the propagation direction is taken into account by incorporating into the detunings the 
Doppler shifts ui c;p v/c, which are practically the same for closely spaced hyperfme compo- 
nents of the ground atòmic state. 

The time evolution of the system's density matrix obeys the equation 

where the operator L accounts for all relaxation processes. In our model they are determined 
by the spontaneous emission from the upper level |3) to the lower levels and the damping 
of lower-level coherence. For a simplicity, all optical decay rates are assumed to be the 
same: 73^=7, and the coherence decay rates of the lower levels are taken equal 7íj— 7o<<7, 
i, j=0,l,2. 

In general case the Eq.(3) cannot be solved exactly. However, an analytical solution 
can be found in symmetric case, when the coupling field is tuned to resonance with |0>^ 
|3>transition in the presence of magnètic field, i.e.A c = 0, while for the probe field the 
condition üü 3í — u p = 0, i— 1,2 holds, so that Ai )2 = —#1,2 = ±5. We consider the case of 
strong coupling làser and weak magnètic fields such that the EIT conditions are satisfied: 

Ql» 5A D , 70(7 + A D ), (4) 

where is the Doppler width of the optical transitions. Then, with the assumptions 
of Q p , 5 << ÍÍ C )7) we are a ble to derive the solution for the Zeeman coherence and the 
absorption a 12 and dispersion /5 lj2 coeficients of probe field in the lowest order of small 
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parameters íl p /íl c and S/íl c . With equal initial populations in the Zeeman sublevels pn(0) 
=P22(0)=0.5 this solution for time t much larger than the upper level inverse decay rate: 
t>>7 -1 , has the form: 



<t) = Y M^), #(t) = ^Re(^E), i = 1,2 (6) 

where 

P3i = + iríO - Ül ^ 2 exp(- 7d í) exp[.( - 6 2 )t] + O(e^) (7) 



P32=P3l(l"2), (8) 



ld = 70 + 2mj57^ (9) 

Here io — — ^> °" — 4:iruj p fi 2 / HcT is resonant absorption cross-section, T = 37/2 is optical 
dephasing linewidth and N is the atòmic density. From the symmetry of the system with 
respect to two circular components of probe light it follows that a%(t) = a 2 (t) = a(t) 
and f3\{t) = — /3 2 (í), as is apparent from Eqs.(6-8). Note also that due to the condition 
Q p << T nearly all of the atòmic population remains in the states |l>and |2>, so that the 
depletion of coupling field can be neglected. 

The first terms in Eqs.(5) and (7,8) represent the results of EIT regime in A- 
subsystems (0-3-1) and (0-3-2) and, of course, coincide with that obtained in phase-sensitive 
model || up to the factor 1/2 emerged from initial conditions. The immediate interaction 
between the two polarization components of probe field in three-level configuration (1-3-2) is 
responsible for the second terms, which reflect the quantum beats in the system. It is worth 
noting that in the limit B— >0 the two ways yield the same contribution to Zeeman coherence 
P2i- We recognize from Eqs.(5-9) that in this limit there is no oscillation in the system. 
At the same time it is seen that the frequency of the oscillations is linear in the magnètic 
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field, while their damping rate 7^ has square dependence on B. As it follows from Eq.(9), for 
2rí2p<5 2 /í2^<<7o the damping of oscillations is determined solely by Zeeman decoherence 
rate 70. In Fig.2 we show the time evolution of the real part of p 2 i obtained by numeri- 
cal integration of Eq.(3) for two different vàlues of magnètic field with the initial conditions 
Pn(0) =p 2 2(0)=0.5 and the other Pij—0 (i,j=0-3). As expected, the oscillations occur in time 
region t>>r _1 . Note that they always start with same value Re(p2i)~ — |^i^2|/(^c + 2^), 
but according to Eq.(9) the amplitude of the oscillations is damped fast as the strength of 
magnètic field increases. In the absence of magnètic field P21 is purely real and reaches very 
quickly the staedy-state value, which is two times larger ( when disregarding the sign) than 
time averaged value of Re(p 2 i) at (see Eq.(5)). It is not difhcult to be convinced, that 
the analytical solution Eq.(5) is represented graphically by the same curves with insignificant 
deviation. For comparison we plot in Fig.2 the Zeeman coherence in the case of A c = — 5q, 
when the symmetry between the two circular components of probe field is violated. How- 
ever, we observe that this has a negligible effect, if magnètic field is small. This analysis 
allows us to conclude that for sub-Gauss magnètic fields our analytical results represent the 
quantitative physics quite well. 

The equations of motion for the propagation of the intensity l p ~ Q p and phases 
$1,2 of the RCP and LCP components of probe field read 

<é + < io > 

where the group velocities v g expressed as 

Vg = c[l + c^-]- 1 ^ AÜ%T- 1 (12) 

are equal for two polarization components. 

The most transparent results are obtained in the case of optically thin media, where 
aL « 1,L is the atòmic celi length. Using the Eqs.(6) and (7) this condition can be written 

as 
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In this case the integration of Eqs.(lO) and (11) yields for the probe transmitted power 
and rotation angle $ = ($ 2 — < ^ ) i)/2 

í%(L,t)=í%(0,t-L/v g ) = íÇ(0) (14) 

L 

/SL íl (0) SL 
P 2 (t - z/v g )dz = — + sin(— ) cos(25í) exp(- ld t) (15) 

Vg i í c f g 







where í~2 (0) corresponds to the input probe field. Here we have used the Eq.(12) for v g and 
have taken into account that 7^ << 5. 

The Eq.(15) shows an interesting variation of oscillating term on parameter r\ = 
SL/vg, which is the ratio of probe propogation time L/v g to the atòmic oscillation char- 
acteristic time S^ 1 . This modulation can be understood rather symply. When rj << 1, 
the àtoms on propogation line are in the same phase of ground-state oscillation during all 
passing of probe field through the medium. This leads to total accumulation of polarization 
rotation and to the maximal relative contribution of oscillating term in output rotation an- 
gle. In opposite case of i] >> 1 the probe pulse moves very slowly and sees the àtoms in 
different phases of oscillation. As a result, the polarization rotation varies periodically with 
distance travelled that eventually leads to the sinusoidal dependence of the total rotation 
angle on the celi length. In this case, of course, the relative contribution of oscillating term 
is reduced, even if sini] = 1. By the same reasons, one can conclude that in high density 
case ( aL » 1) the oscillating term in <E>(i) is suppressed even more. 

The estimations with the help of Eq.(15) show that in the case of D 2 line in Rb 
the polarization rotation as high as 0.1 rad can be achieved for sub-Gauss magnètic fields 
(B~0.7G), atòmic density N~10 10 cm~ 3 , L=3cm and íl c = T (aL~0.1) with the amplitude of 
oscillations ííp/ííj?~0.1. For large vàlues of rj, we should make sure of compatibility of the 
conditions i] » 1 and aL << 1. Combining the parameters a,i] and using the expression 
of v g from Eq.(12) the latter condition can be written in the form 2f]ST/íll << 1, which can 



be fulfilled in a wide range of 77, if one takes into account that 5 « í7 c .The predicted effect 
is much easy to be observed in a cooled atòmic ensemble, where the Doppler broadening can 
be neglected that allows to avoid the strict restrictions of Eq.(4) on the magnètic field. 

It is worth noting that the present technique can be used to fully entangle a pair of 
very weak optical fields. There have been several studies |TB| . |T4| predicting an unprecedental 



strength of nonlinear interaction between two low-power light pulses. The attempts to 
apply these schemes to generation of photon entangled states [Ï5| have confronted with the 



problem of the mismatch between the slow group velocities of propagating weak fields that 



restricts their effective interaction time ||16|| . To overcome this difficulty it was proposed 
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Ï5| to use the mixture of different isotopes of alkali àtoms and to apply an appropriate 
magnètic field. To demonstrate how the efficient interaction of two weak pulses arises in 
our scheme we consider the case of two uncorrelated weak beams with equal freqüències 
instead of two circular components of linearly polarized probe light and suppose that the 
frequency of the fields is chosen such that, say, the first field is in exact resonance with the 
|1>— ► |3>transition in the presence of magnètic field, i.e. Aj = 0, while A 2 = 25. Then, 
when putting in Eq.(7) 8\ = and 5 2 = 25, we find that the phase shift and absorption 
of the first field are proportional to the intensity of the second one that creates the cross- 
modulation between weak beams. The main advantage of our scheme is that the problem of 
group velocity mathching does not arise here, since, as it was shown above, the two beams 
have the same group velocities and, hence, the interaction between them is maintained for 
a very long time. 

In summary, we have described new phenomenon of long-lived temporal oscillations 
of NFR, which arise when the conventional Faraday configuration is driven by a strong 
làser field, and have shown that these oscillations reflect the quantum beats in coherent 
superposition of ground-state Zeeman sublevels. We have found the analytical solution for 
total angle of Faraday rotation,which allows to understand quantitatively the main features 
of the effect. Our study predicts also the ability of the proposed scheme to achieve an efficient 
creation of two-photon entanglement. This question will be discussed in detail elsewhere. 
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Figure captions 

FIG.l. Four-level model for observation of oscillating Faraday rotation. íl c and íli ; 2 are 
the Rabi freqüències of coupling field and a± components of probe light, respectively. A 2; i 
are the shifts of m=±l levels induced by magnètic field. 

FIG.2. Numerical results for time evolution of Zeeman coherence Re(p 2 i) i n the case of 
Q C =T, ííi )2 =0.3r, 7o=10" 4 r, A c = and for different vàlues of magnètic field B: 1) 5 = 0, 2) 
5=0. 03r and 2) S—0. 1T. The dotted line corresponds to the case of A c = 5=0. 03T with the 
same vàlues of the rest parameters. The Doppler broadening is ignored. 
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FIG.2 




